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Motivated by recent attempts to reconcile hints of direct dark matter de- 
tection by the CoGeNT and DAMA experiments, we construct simple par- 
ticle physics models that can accommodate the constraints. We point out 
(— I \ challenges for building reasonable models and identify the most promising 

Qh' scenarios for getting isospin violation and inelasticity, as indicated by some 

pi • phenomenological studies. If inelastic scattering is demanded, we need two 

O I new light gauge bosons, one of which kinetically mixes with the standard 

model hypercharge and has mass < 2 GeV, and another which couples to 
baryon number and has mass 6.8 it GeV. Their interference gives the de- 
sired amount of isospin violation. The dark matter is nearly Dirac, but with 
small Majorana masses induced by spontaneous symmetry breaking, so that 



0^ ■ the gauge boson couplings become exactly off-diagonal in the mass basis, and 

the small mass splitting needed for inelasticity is simultaneously produced. 



If only elastic scattering is demanded, then an alternative model, with inter- 
' ference between the kinetically mixed gauge boson and a hidden sector scalar 

Higgs, is adequate to give the required isospin violation. In both cases, the 
light kinetically mixed gauge boson is in the range of interest for currently 
running fixed target experiments. 

X' 

: 1. Introduction 

Hints of direct detection of dark matter (DM) currently exist from two experiments. 
There is a long-standing observation of an annual modulation in the signal observed by 
DAMA [1], whose statistical significance is beyond question. Last year the CoGeNT 
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experiment reported excess events in their lowest electron energy bins [2] , followed more 
recently by a 2. 80" detection of annual modulation in the signal [3]. Under the simplest 
assumptions about the nature of the dark matter interactions, these two observations 
appear to be incompatible with each other [U [5] and with upper limits obtained by 
other experiments, especially CDMS [HIT], XenonlO [8] and XenonlOO [9j. Channeling 
of recoiling ions along the crystal planes in the detectors has been suggested as one 
loophole for reconciling the conflicts, but this has been argued to be too small an effect 
by ref. [TU]. Uncertainties in quenching factors can also be used to help reconcile the 
two positive detections [TT]. 

There are alternative microphysical ways of ameliorating the tensions that have been 
explored in recent papers [I2]-[T9]. One important modification is to allow for isospin- 
violating interactions of the dark matter with nucleons, since the coherent spin-independent 
matrix element is a sum over the interactions with protons and neutrons: 

aN-^{fpZ + fn{A-Z)fan (1) 

where A, Z are the number of nucleons and protons in the nucleus, respectively, and 
an,N are the respective cross sections for DM scattering on neutrons and the nucleus. If 
isospin is conserved then fp = fn, but if fn/ fp ~ —0.7, the limits placed on an by the 
Xenon experiments are greatly relaxed \12\ IT3j . 

Moreover, several groups have indicated that inelastic scattering of DM with a 
small mass splitting (5 ~ 10 keV can have a beneficial effect either for the the agreement 
between DAMA and CoGeNT [l5], the conflict between CoGeNT and CDM^p2lfT6]. 
the goodness of fit to the CoGeNT modulated signal alone [18j, or marginally the overall 
goodness of fit to all data p9j . 

In this work, we look for the simplest and most natural kind of models that could ac- 
commodate these two generalizations, isospin violation and inelasticity, with an emphasis 
on hidden sector models with Higgs or gauge kinetic portals to the standard model. (For 
previous model building efforts which do not focus on these aspects, see references 12U]- 
[28].) A similar study to the present one was done in ref. [29], but considering more 
elaborate models than we examine here. Our inelastic model also has elements in com- 
mon with that of ref. |30j . although the latter incorporated neither isospin violation nor 
inelasticity. (See [31] for another recent isospin-violating model.) In ref. [32], isospin 
violating couplings of a single vector were obtained through a combination of kinetic 
mixing and mass mixing with the Z'. We will argue that interference between two new 
vector mediators is an elegant way of getting isospin violation if one demands that only 
inelastic scattering takes place. It is well known that Dirac states coupling to vectors 
become off-diagonal in the mass eigenstate basis if small Majorana masses are introduced 
|14j . Thus it is not challenging to account for the inelastic nature of the couplings. 

On the other hand, ref. [TS] finds that inelasticity improves the global fit to the data 
only moderately, so that one might also contemplate models with no DM mass splitting 

Inelastic scattering increases the ratio of the modulation amplitude to the unmodulated rate, softening 
the discrepancy with CDMS W, which reports no evidence of events that would be compatible with 
the CoGeNT signals. 
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and purely elastic scattering. The overall fit is in fact not very good, indicating either 
that some of the data have inconsistencies or that the dark matter interpretation is not 
correct. In this work we will assume the former, in which case one might be motivated 
to consider inelasticity as a secondary criterion, which might or might not survive as 
the data improve. Accordingly, we consider both elastic and inelastic models in the 
following and leave it to the reader to judge how strongly the latter is preferred by the 
data. Elastic scattering naturally arises in alternative models where the isospin violation 
comes about by interference between vector and scalar exchange. We will construct our 
models for isospin-violating dark matter starting with this simpler possibility, showing 
why it does not naturally accommodate inelastic couplings. 

2. Scalar versus vector exchange 

A natural way to induce scalar-mediated interactions between the dark sector and the 
standard model (SM) is to introduce a new Higgs field (p that is a singlet under the 
SM gauge group, and communicates to the SM through the interaction XcfP'h?'. If (f) 
gets a VEV then the mass eigenstates are admixtures of (p and h with mixing angle 
9s = A(^)u(m^^ — 771^^), where m/i^^ stand for the mass eigenvalues (we assume small 
mixing, 0^ <C 1) and v = 246 is the SM Higgs VEV. If the DM x has a Yukawa coupling 
yx4'X^ then cj) mediates interactions with SM fermions / with the strength yy fd s / m^, 
assuming that <^ nifi. Here y/ is the SM Yukawa coupling of the Higgs to fermion 
/• 

In this scenario, cj) couples to SM matter proportionally to h. Although these couplings 
violate isospin, it is a very small violation when it comes to scattering on nuclei, because 
nucleons get almost none of their mass from the valence quarks. The couplings of the 
Higgs to nucleons are dominated by the sea-quark and gluon content, which are the 
same for neutrons and protons (see for example j23|l29j). Therefore we get no significant 
isospin violation from couplings mediated by the Higgs portal: it has fp = fn- 

If in addition we invoke a gauge kinetic mixing portal, {e/ cos 6w)F^i,Z'^'^ (where F 
and Z' are the respective field strengths of the SM hypercharge and the hidden sector 
U(l)', and 9w is the Weinberg angle), then interference between the new vector and the 
new Higgs does lead to tunable isospin violation. This is because the vector only mixes 
significantly with the photon, leading to fn/ fp = 0. Of course, this is not the ratio we 
need for DAM A and CoGeNT. But interference between the vector and scalar exchange 
allows for any desired value. If g' is the new U(l)' gauge coupling to the DM, then 

f_n ^ {yyn9slm\) 

fp {g'ee/ml,) + {yyn9s/ml) 

where y„ is the coupling of h to the nucleon: y„ = 0.36 m„/f in terms of the nucleon 
mass rUn and the Higgs VEV v [33l |3l] . 

The main objection to this scenario arises if we want the couplings to DM to also be 
inelastic. Suppose there are two mass eigenstates x± with masses M± = M it (5/2 split 
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by the small amount 6. Let us first write down an effective potential: 

y = ^ X±M±x± + \^m\,Z''^ + + gx+$' X- + yx+(t>X- (3) 

(It is also understood that Z' couples to the electromagnetic current with strength ee and 
(j) to the SM fermions / with strength OsVf-) The mass splitting can arise by spontaneous 
symmetry breaking through couplings of the form (j)*{yLXLXL + yRXRXR) where Xl,r are 
the Weyl components in the original Lagrangian, and (j) carries twice the U(l)' charge 
of X- 

The troublesome question is how the scalar interactions came to be purely off-diagonal 
in the mass eigenbasis, given that (j) has a VEV which is needed in order to mix with 
the SM Higgs h. To arrive at ([3]), there must have been a bare mass term fix+X- (here 
expressed in the mass eigenbasis) that was exactly canceled by Iii the ab- 

sence of this tuning, (j) will have diagonal plus off-diagonal couplings, and this presents a 
complication for getting the desired level of isospin violation, since the interference be- 
tween (j) and Z' only occurs in the inelastic channel. This is because the Majorana vector 
couplings are purely off-diagonal (see the following subsection) . We do not contemplate 
this finely-tuned situation any further. 

However if instead of couplings of the form </'*(xi,Xi, + XrXr) ™- which </> must carry 
a compensating charge, we have the interactions 

X (M + y4>) X (4) 

in terms of a Dirac fermion x, then x remains Dirac after (j) gets its VEV, and there is 
no mass splitting, (f) can interfere with Z' through the purely elastic couplings. Even 
though we must invoke an additional singlet Higgs (j) to break the U(l)' symmetry, since 
now (j) is neutral under U(l)', this is still an economical model, whose consequences we 
will consider. For simplicity we take (j) to be real. We do not consider the case of a 
pseudoscalar (275) coupling because this leads to a nuclear scattering amplitude that 
is suppressed by the DM velocity, making it more difficult to interfere with the vector 
exchange contribution to get the desired isospin violation (since the latter has no such 
velocity suppression). 

2.1. Nondiagonal gauge couplings 

In contrast to the couplings of a scalar to DM, the purely off-diagonal coupling of the 
gauge field is natural if x± are Majorana fermions that originated from a Dirac particle 
before spontaneous symmetry breaking [14j. Consider the interactions (again in the 
model where (p carries two units of the x charge) 

V = \xlMxr + {Xl Pl Xl + Xr Pr Xr) + h.c. (5) 

where now = i''pL,<^2''Pt)^ X^ — {~^2'^%'4'r) denote Majorana-Dirac spinors con- 
structed from the Weyl spinors, here renamed ipL,R to avoid confusion. When (j) gets a 
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VEV, the mass matrix becomes 



( 



M 
M n 



) 



(6) 



which is diagonaUzed by Xl,r = ~^(.X+ =t X-) with mass eigenvalues \M±\ = M ± fi, 
where fi = y{(j)) <^ M. If the U(l)' interaction was orginahy vector-hke, then it becomes 
exactly off-diagonal in the mass basis because there is no vector current for a single 
Majorana state. This is a strong motivation to prefer vector mediators if we aim for 
both isospin violation and inelastic off-diagonal couplings. 

3. Interfering vector exchanges 

The previous discussion motivates us to build a model in which the interfering scalar 
current is replaced by another vector current. The new vector need only couple to 
isospin differently from the kinetically mixed U(l) that has fn/fp = 0. Coupling to 
B — L is attractive from the point of view of anomaly cancellation, but such couplings 
are very strongly constrained because of the leptonic interactions (see for example [35j). 
The simplest possibility that avoids these constraints is coupling to B alone. U(1)b 
is anomalous and it also has mixed anomalies with the SM gauge groups, that can 
be canceled by adding the appropriate exotic heavy particles [36] -[II]- We will not 
discuss the implications of these new particles further here, although they can provide 
complementary collider signatures to test the model. Our addition of a single vector-like 
DM particle coupling to B does not spoil the anomaly cancellation achieved in these 
models. 

We refer to the \]{1)b gauge boson as B^ and for simplicity assume that it couples 
with equal strength qb to the DM and to the SM baryons. It also couples with equal 
strength to protons and neutrons, just like the singlet Higgs of the previous section. 
Therefore it is clear that ([2]) is replaced by 



This model is almost complete, but we have accounted for the breaking of only one 
linear combination of the two new U(l)s through the VEV of (p. Notice that (f) must 
have charges —2{g',gB) under 11(1)' x U{1)b in order for ([5]) to be gauge invariant. To 
completely break the symmetry we need another field <j) with different charges. This 
means that the mass eigenstates for the gauge bosons are generally admixtures of the 
original fields, and that both will therefore kinetically mix with the SM hypercharge. 
We need to clarify the relation between the couplings appearing in ([7]) and the original 
Lagrangian parameters. 

A simple way to ensure that the above relations are approximately correct despite 
mixing of the new gauge bosons is to assume that 3> by assigning cj) the 

charges (O,^^) such that ciBii')'^ ^ 9b {4')'^ ■ In that case the mixing is suppressed by the 



fn ^ {gl/ml) 



(7) 



fp (g'ee/m|,) + {gl/m%) 
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Model 


cr„ (cm^) 


M (GeV) 


6 (keV) 


fp/ fn 


vector Bfj^ exchange 


3 X 10"^^ 


8 


9.3 


-1.53 


scalar (p exchange 


6 X 10-39 


7.5 





-1.54 



Table 1: Best-fit values of the DM-neutron elastic scattering cross section, DM mass, 
mass splitting, and isospin violation, from ref. [19], appropriate to the given 
theoretical model. 

large mass differ enc^ and we can take 

mB = ~gB{4>) > mz'^2g'{cl)) (8) 

The mixing angle is approximately 9 = {gs / g'){fnz' / ''t^b)'^ ■ We will find that its value 
scales proportionally to the gauge kinetic mixing parameter e, such that = Ae. (This 
relation follows from eq. ()10p below and the relic density constraint, fig. [lU^a).) 

4. Determining the couplings 

We must show that values for the parameters exist that can give the right cross section for 
CoGeNT and DAMA, and the right relic density for the dark mattero For the effective 
elastic cross section of DM on the neutron, the DM mass and mass splitting, and level of 
isospin violation, we will consider the three cases shown in tabled! which is corresponds 
to the best-fit values found by ref. [TO] for the cases of endothermic X-^ ^ i aiid 
elastic scatterings, respectively. These are the ones appropriate to our two models. We 
will show that the exothermic reactions, ~^ X-^ ■• ^t-re possible when e ^ IQ-^'^, 

whereas otherwise the excited state is depleted by x+X+ ~^ X-X- downscatterings in 
the early universe. 

The theoretical cross sections for DM-neutron scattering in our models, in the elastic 
limit, are 

o- = ^ X I 9%/^%^ exchange | 

vr \ [yUnOsf /m^^-, scalar (/> exchange J 

where is the reduced mass for the DM-nucleon system and v is the DM velocity. By 
equating the (T„ values in tabled] to those in eq. ([9]) and using eqs. ()2l7p . we obtain 

2 

— = 232 GeV, ^ = -(79.9 GeV)^ (B^ exchange) (10) 
9B 9e 

2 

- 346 GeV, ^ = -(118.9 GeV)^ (scalar (j) exchange) (11) 



^If {(j)) — u and {(p) = then the gauge boson mass matrix in the basis (B, Z') is ^ ^^4g '^g^S^ ^Ig'^-u^ ) 
^ Since our DM is Dirac or quasi-Dirac, there is the interesting possibility for an asymmetry between x 
and X being the origin of the rehc density, which has been widely discussed in the recent Uterature 
(see for example ^42,). For this work we will assume the asymmetry vanishes. 
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To get the correct mass spectrum for the model with i?^, we can set the bare Dirac 
mass directly to M = 8 GeV, and choose y{(j)) = 4.7 keV. Recall that only Z' gets its 
mass primarily from 0, in this model, so if e is sufficiently small, it is possible to have 
{(p) ~ 10 GeV or less. The Yukawa coupling still needs to be quite small in that case, 
y ~ 0.5 X 10~^. However this is only 4 times smaller than the electron Yukawa coupling 
in the standard model, so it is not unreasonable. Alternatively, for the purely elastic 
models with cj) exchange, we need to set M + y{(j)) = 7.5 GeV. 

4.1. Relic density 

Next we consider the relic density. Starting with the vector exchange model, there 
are several possible annihilations into gauge bosons: X±X± ~^ BB^ X±X± ~^ Z' Z' , 
X±X± ~^ BZ' . In addition there are coannihilations into quarks, x+X- ~^ QQ mediated 
by B in the s-channel, and also x+X- ~^ ff ™to all charged SM fermions / except the 
kinematically inaccessible top, mediated by the Z'. Averaging over all the possibilities, 
we find the annihilation cross section 

(o-ann-f) = j^'^ S {g' , gs , e, x' , xb) (12) 

where 5 is a dimensionless function of the couplings and the mass ratios x' = niz'/M, 
xb = itt-b/M, given by 

S = \g'^ h{x' ,x) + g^g'"^ f2{x ,xb) + \g% h{xB,XB) 

+ g%Y. ^c,^h{xB,Xi) + {g'eef N,,iQU^{x\xi) (13) 

i=u,d,s,c,b i= , 

u,d,s,c,b 

Here Xi = rrii/M for SM fermion i, with charge Qi and number of colors N^^i, and the 
functions fi are given in the appendix. 

For the the scalar exchange model, all the analogous processes to the previous case 
are present, with replaced by (j). The XX ~^ 4>4> contribution is p-wave suppressed, 
as is XX ~^ ff exchange, so we neglect them. The xX ~^ ^'Z' and Z'-mediated 
XX ~^ ff contributions are the same as in (|13p . We find that (jlSp is replaced by 

S^^g'^ fi{x', x') + y^ g'^ /4(x', x^) + {g'eef ^ N,,i QI /3(x', x,) (14) 

where x^ = rrifff/M and /4 is defined in the appendix. 

For each of the models, we equate (jl2p to the standard value of the cross section for 
the observed relic density, {(Ta,nnv) = 3 x 10^^^ cm^/s. This gives the constraint 

5 = 1.7x10^5 (15) 

To determine the parameters satisfying (jlSp in the vector exchange model, we assume 
several choices for e that can be compatible with laboratory bounds on kinetic mixing 
(see next section), and use (jlOp to eliminate g',gB in favor of mz',mB- For the scalar 
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Figure 1: (a): Contours in the rriB-mz' plane that give the observed relic density in the 
vector exchange model, for gauge kinetic mixing parameter e = 10~^, 10~^, 
and 10~^. The ellipse highlights the region where the gauge boson mixing 
angle 9 is small as needed for consistency. 

(b)-(d): Analogous contours in the m^-niz' plane for the scalar exchange 
model, for e = 10~^ (b), 10^'^ (c), 10~^ (d), and several values of the Higgs 
mixing angle 9s- 



exchange model, we similarly use (jlip to eliminate y in favor of m^. This case has the 
free parameter of the Higgs mixing angle 9s to be varied in addition to e. For m,^ < 5 
GeV 9s is constrained to be less than 0.01 from the width of the Z boson due to decays 
Z — )• (j)ff and from B meson decays B — )• (pff- |43j . 

We thus obtain, for each value of e, the contour in the uiz'-mB plane corresponding 
to the observed relic density for the vector exchange model, shown in fig. [T][a). Similarly 
in the scalar model, for each pair {e, 9s}, we find a contour in the mz'-rriff, plane, shown 
in fig. [I]I|b)-(d). It is clear from fig. [T]that mz' and tend to be < 1 GeV and both 
scale as y/e. Only the vector B^ can remain somewhat heavier as we now explain. 

Recall that we previously made a simplifying technical assumption, ^> mz', to 
ensure small mixing 9 between the two gauge bosons. This assumption might be relaxed 
somewhat, but at the risk of increasing the highly constrained couplings of B to leptons 
due to the gauge kinetic mixing. To the extent that 9 is small, only the vertical part 
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of the contours where uib = 6.8 GeV is relevant, giving a sharp prediction for the 
mass, if M is known. Since the determination of M could well be uncertain by ±1 GeV, 
we find an uncertainty of ioi GeV in by varying M. In this vertical branch, the 
XX annihilation cross section is dominated by the gj^ contributions in ()13p . 

4.2. Relative abundance of excited state 

In the vector exchange model where we have a small DM mass splitting, the process 
X+X+ ~^ X-X- mediated by the Z' (and the B, although we find the former dominates) 
efficiently depletes the x+ population in the early universe over part of the allowed 
parameter space. The downscattering cross section for a similar model was calculated 
in [2^, which adapts to the present case as 

^ ( g'^ g| y f25_ 

~ \ml, + 2M5^ ml + 2M5) An S M ^ ^ 

These interactions freeze out at temperature Tf given by n{(T^v) = H where n ~ (7 x 
10~^°GeV/Af)r| is the DM number density and H ~ T^Mp is the Hubble constant. 
If the DM remained in kinetic equilibrium with the SM down to Tj then the relative 
abundance of x+ to X- would be suppressed by ~ exp(— J/Ty). However it is the kinetic 
temperature of the DM which is important here, and if kinetic decoupling occurs at a 
temperature Td > Tj, then the suppression is more severe, ~ exp(— (5Td/T|) (see [26] for 
a discussion of this issue.) The kinetic equilibrium is controlled by the scattering of DM 
on electrons, whose cross section is approximately 




M= (17) 



We find that this goes out of equilibrium at Ta = 10^5 = 10 MeV. 

Using this methodology, we estimate that the relative abundance of x+ is unsuppressed 
for e > 10-3-^ The suppression turns on exponentially fast as a function of e, with 
e = 10~^ giving a relative abundance of ~ exp(— 10"^), while e = 10~^ leads to almost 
no suppression, e^*^'^. Therefore it is possible to realize the exothermic dark matter 
scenario suggested in ref. over some part of the allowed parameter space. In these 
cases the scatterings will be an average over the endothermic and exothermic ones since 
both states are equally populated. We note that ref. [19j finds a moderate preference for 
exothermic reactions in their global fits. 



5. Discussion 



To recapitulate, we have investigated two hidden sector dark matter models that violate 
isospin in an optimal manner for reconciling the CoGeNT and DAMA signals with 
constraints from XenonlO and XenonlOO. The vector exchange model has two Majorana 
mass eigenstates x± (or ^ single pseudo-Dirac particle) with masses M± = M ziz 5 = 8 



9 



GeV ±5 keV, and two light gauge bosons Z' and B with masses m^/ < 2 GeV, ms — 6.8 
GeV, and the couplingqfl 

X+{g'^' + 9B^)X~ + gnfeB^ + ef^AiK + ^^m) (18) 

Here e < 10~^ is the gauge kinetic mixing parameter for the Z', ^? = 4e is the mixing 
angle between B and Z', and j^^,^^ are the respective electromagnetic and baryon 
number currents of the standard model. The couplings are adjusted to give the optimal 
isospin violation fp/ fn — —1.5 [19j via eq. ([7]). There must also be diagonal Yukawa 
interactions yx±4'X± to a singlet Higgs (j) whose VEV leads to the small mass splitting, 
but y ~ 10~^ is much smaller than the gauge couplings (qb = 0.029) and therefore 
we have neglected these interactions. The Dirac mass M = 8 GeV appearing in the 
original Lagrangian is protected by chiral symmetry, and so does not introduce any new 
hierarchy problem. The small scale ((/>) ~ GeV on the other hand is unexplained unless 
one invokes supersymmetry in the hidden sector [IS] or some other UV completion. 

The scalar exchange model has exactly Dirac dark matter with M = 7.5 GeV (hence 
only elastic scattering) and a real singlet (J) that mixes with the SM through the Higgs 
portal XcfP'h'^. Like the previous model these also have the kinetically mixed Z' vector. 
The interactions are given by 

x{g'i' + #)x + esY,y^Uf^ + ^JekK (19) 

i 

where fi are the SM fermions with their Yukawa couplings We find that ^ 10 
GeV to satisfy the relic density constraint. If happens to be close to this upper limit, 
it could be discoverable at the LHC, while for the very light cases < 1 GeV, indirect 
discovery could come from rare decays such as -B — )• (j)X followed by — )• [43]. 

One of the most exciting aspects of these models is that they predict new low-energy 
interactions mediated by the light Z' with a strength relevant for detection in beam- 
dump experiments [461 HZj such as APEX [48J and the Mainz Microtron [39| and the 
low-energy e~^-e~ collider experiment KLOE |50j. The still-open window of parameter 
space in the e-mz' plane corresponds roughly to that which we have identified in this 
paper as being compatible with the x± relic density. In our model, m^/ is only bounded 
from above, depending on the value of e, as shown in fig. (T) m^/ <, 2y^e/lCp^ GeV in 
the vector exchange model. 

On the other hand, the gauge boson B of baryon number is extremely hard to detect 
due to its very weak coupling and relatively large mass. For example, constraints from 
new contributions to T decay into quarks are easily satisfied [51]. The Tevatron sets 
limits on gs ^ 0.6 from the nonobservation of pp — ?■ — )■ XXJ where j is a single jet 
[52j . which is also satisfied by our model. Ref. [37| pointed out that the kinetic mixing of 
B leads to weak Tevatron constraints from the Drell-Yan production of lepton pairs. It 

^There are also couplings to the Z boson current given by — e tan 9w JTiJ^ jz ("^z' Z'^ + 6 m% B^) 
that come from the mixing of Z' with weak hypercharge; see for example 144:. Larger contributions 
to the kinetic mixing of the B with Z can be generated from SM loops below the scale of baryon 
symmetry breaking [37) . 
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may be interesting to update these constraints since [37] was written before the upgrade 
of the Tevatron. 

The best indirect confirmation of its presence will be the discovery of an exotic extra 
family of quarks with baryon number ±1 [SQJ- In the simplest such models \37\ I38j. this 
fourth generation gets its mass through the usual couplings to the Higgs, requiring its 
mass to be at the electroweak scale and limited by large Yukawa couplings leading to a 
Landau pole near the TeV scale. But ref. [39] shows that this limitation can be removed 
using vector-like quarks (from the point of view of the SM SU(2) gauge symmetry) and 
giving mass to them through the VEV of the field which breaks baryon number, cp in 
our model. If the coupling gs in ^ is sufficiently small, for example gB ~ 0.007, then 
((/)) can be at the TeV scale. 

Astrophysical constraints are rapidly closing in on light dark matter models. If ^ 10 
GeV DM annihilates predominantly into e^'^e^ with the standard relic density cross 
section, it is ruled out by its distortions of the CMB [531 El]- The n'^iJ,~ channel is 
still open since a large fraction of the muon energy is converted to neutrinos which have 
no effect on the CMB. Other channels have an intermediate effect between these two 
extremes [53J. In our model with two vector bosons, the annihilation is primarily into S's 
followed by decay into light quarks, which appear to be still be allowed. On the other 
hand, models that produce too many neutrinos are constrained by SuperKamiokande 
limits on XX ~^ from the sun |55j . Even in optimally isospin violating models such as 
we have considered here, annihilation of ~ 10 GeV DM with the required cross section on 
nucleons for CoGeNT/DAMA is ruled out for the bb channel and marginally allowed for 
cc and lighter quarks as in our vector model. These limits can be improved in the future 
using data from IceCube/DeepCore [56]. Finally, annihilations of light DM in dwarf 
satellite galaxies of the Milky Way that produce too many gamma rays in association 
with charged particles have recently been severely constrained by Fermi data |571 [58] 
(see also [59]). Again, the bb channel is excluded but annihilation into light quarks is 
still allowed. 

After the first version of this paper was posted, ref. [6^ appeared, which considers a 
similar class of models. 

Note added: After completing this work we became aware of refs. |6H [62] showing 
that PAMELA antiproton constraints are in conflict with a B vector boson of mass 
greater than 2mp, which would favor the lower-mB parts of the contours of fig. 1(a). 
In addition, we discovered that the gauge boson mixing angle effects in the present 
model cannot be ignored even when 6 is small, due to the occurrence of l/m\i in the 
amplitudes, which scales as 1/9. This can be overcome by introducing an additional 
contribution to ms- We intend to address these issues in a forthcoming publication. 
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A. Kinematic functions for annihilation cross section 

The functions of mass ratios appearing in the annihilation cross sections are as follows. 
They were computed using Feyncalc [63\ . 

MXI,X2) = ) — ^e(l-X2), 

(1 - h^l) 



f2{xuX2) = ^ ^^^^ / ' 2 e(2-Xi-X2), 



'^-\{xl + xl)Y 
/3(a;i,X2) = ^ , 1 9x2 6(1-^:2) 

(1 - -A) 



fi{xi,X2) 



\ + x\ ix2 + Tr{x\ 



X: 



2\2^ 



[1 \{x\ + X^) + j^{x\ X^)^) 



2\2\l/2 



'l-\{xl + xl)Y 



e(2-xi -X2) 
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